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In this presentation I summarize recent work on Cepheid research carried out at 
the National Central University. The mid-infrared period-luminosity (P-L) relations for 
Cepheids are important in the James Webb Space Telescope era for distance scale work, 
P> . as the relations have potential to derive the Hubble constant within ~ 2% accuracy - 

a critical constraint in the precision cosmology. Consequently, we have derived the 
mid-infrared P-L relations for Cepheids in the Large and Small Magellanic Clouds, 

■ using archival data from the Spitzer Space Telescope. Kepler Space Telescope is a 
NASA mission to search for Earth-size and larger planets around Sun-like stars, by 

Q , observing continuously the stars in a dedicated patch of the sky. As a result, the almost 

^ • un-interrupted observation is also used for stellar variability and asteroseismological 

' study. However, Kepler observations are carried out with a single broad-band filter, 

. hence ground-based follow-up observation needed to complement Kepler light curves 

to fully characterize the properties of the target stars. Here I present the ground-based 
optical follow-up observations for two Cepheid candidates located within the Kepler's 
field-of-view. Together with Kepler light curves, our ground-based data rule out V2279 

■ Cyg being a Cepheid. Hence VI 154 Cyg is the only Cepheid in the Kepler's field. 
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1. Introduction 

Classical Cepheid variables (hereafter Cepheids) are post main-sequence yellow su- 
>v>( \ pergiants. Cepheid masses range from ~ SM© to ~ 11 M© with surface temperatures 

j_j ■ between ~ 5000^ to ~ 6500^. Cepheids obey the period-mean density relation, hence 

5^ \ their pulsating periods (1 < P < 100, where P is period in days) are well corre- 

lated with luminosity. This is known as the period-luminosity (P-L) relation, also re- 
ferred as the Leavitt Law, commonly written in the form of = ax log(/') -i- b;^ or 
mx = ax log(P) + bx- Since Cepheids are evolved pulsating stars, they are ideal labora- 
tories for testing the stellar pulsation and evolution theories. The Cepheid P-L relation 
is an important rung in the cosmic distance scale ladder that can be used to calibrate 
a host of secondary distance indicators (for examples, the Tully-Fisher relation, peak 
brightness of Type la supernovae and surface brightness fluctuations). These secondary 
distance indicators, located well within the Hubble-flow, can then be used to measure 
one of the most important parameter in modem cosmology - the Hubble constant. 

In this proceeding, I will present some recent progress on Cepheid research car- 
ried out at the National Central University (NCU), that focused on the following two 
topics: (a) derivation of the mid-infrared Cepheid P-L relations based on archival data 
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from Spitzer observations; and (b) investigation of potential Cepheids located within 
the Kepler field using ground-based follow-up data. 



Derivation of the Mid-Infrared Cepheid P-L Relation Using Spitzer Archival 
Data 



In the era of precision cosmology, it is desirable to independently derive an accurate 
Hubble constant, This is because degeneracy exists in CMB (cosmic microwave back- 
ground) anisotropics between_the flatness of the Universe and the Hubble constant. 
Therefore, as stated in lHul (l2005h : 



"... to test the cosmological constant hypothesis and measure the equation 
of state of the dark energy at z ~ 0.4 - 0.5, the best complement to current 
and future CMB measurements is a measurement of the Hubble constant 
that is accurate at the few percent level." 



Furthermore, Figures 23 and 24 from lMacri et al.l ( 20061) clearly illustrate the improved 
precision of measured dark energy parameters when the accuracy of Hubble constant 
is increased. Current systematic errors of Hubble constant measurements based on 
optical data, using the Cepheid P-L relation as first rung of the distance scale ladder, 
is ~ 5% to ~ 10% (see the review given in iFreedman & Madorell20l"ol) . In the near 
future, it is possible to reduce the systematic error of Hubble constant to ~ 2% from 
mid-infrared obser vations from the Spitzer an d/or next generation James Webb Space 
Telescope (JWST) dFreedman & Mador3l2010), taking huge advantage of the fact that 
extinction is negligible i n the mid-infrare d /preedma n et all l2008l : Ngeow & Kanbud 
12008 : Ngeow et al..2009l : IFreedman & Madore 2010). The first step toward this goal is 
to derive and calibrate the Cepheid P-L relation in the mid-infrared. 

2.1. The IRAC Band P-L Relations 

Prior to 2008, the longest available wavelength for the Cepheid P-L relation is in the 
K band. In 2007, archival data from SAGE (Surveying the Agents of a Galaxy's Evo- 
lution) was released, where SAGE is one of the Spitzer Legacy programs that map 
out the Large Magellan ic Cloud (LMC) in Spitzer IRAC bandj] (and MIPS bands, 
see iMeixner et al.l 2006h . Hence, mid-infrared Cepheid P-L relations can be derived 
by matching the known L MC Cepheids to the SAGE catalog - this has r esulted two 
papers pubhshed in 2008 (IFreedman et all l2008l : iNgeow & Kanbutll2008h . Both pa- 
pers utilized the SAGE Epoch 1 data matched to different samples of LMC Cep heids: 
^eedman et al. (2008) matched to ~ 70 Cepheids from iPersson et al.l (Hool) while 
Ngeow & K anbur (2008) used the OGLE-II Catalog (Optical Gravitational Lensing Ex- 
periment,[Udalski et al. 1999) that includes ~ 600 LMC Cepheids. 

After the release of SA GE Epoch 1 an d 2 d ata, the IRAC band L MC P-L relations 
have been updated in Madore et al.' ('2009') and 'Ngeo w et all (l2009i) by averaging the 
two epochs photometry. Ngeow et al. (2009) also included the ~ 1800 LMC Cepheids 
from the latest OGLE-III Catalog ( Soszynski et al.ll2008l with ~ 3x more Cepheids 
than OGLE-II Catalog). In addition to LMC, IRAC band P-L relations have also been 



'The IRAC bands include S.d/im, 4.5yum, 5.8;um and S.Oyum bands. 
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derived for SMC Cepheids (see iNgeow & Kanbid 120101 . for details) using the pub- 
licly available SAGE-SMC data. Finally, IRAC band P-L relations for a small number 
of Galacti c Cepheids that possess independent distance measurements have been pre- 
sented by iMarengo et alj ( 2010h . All together, there are six sets of IRAC band P-L 
relations for Cepheids in our Galaxy and in Magellanic Clouds. Slopes of these P-L 
relations are summarized in Figure [T] 



□ GALl ■ LMCl 
" GALS ■ LMCa 
o GAL3 • SMC 



a.y 8.4 

ia+iog(o/H) 




ia-iog(o/H) 



Figure 1. Slopes of the empirical IRAC band P-L Relations from three galaxies. 
GALl and GAL2 were derived from Galactic Cepheids with infrared surface bright- 
ness distances using the "old" and "new" p-factor, respectively; while GAL3 was 
based on 8 Galactic Cepheids with HST parallax measurements (for more details, 
see Marengo et al. 20 10). LMCl a nd LMC2 ai-e the empirical LMC P-L slopes from 



iMadore et all (I2009I) and iNgeow e t al. (2009), respectively. SMC P-L slopes are 
adopted from lNgeow & Kanbud (12010) . These P-L slopes can be grouped into the 
"shallow" slopes and "steep" slopes. The Horizontal dashed lines represent the aver- 
aged slopes in these two groups, with the values given on the left of these lines. The 
dashed and dotted lines are the Icr boundary of the averaged values for the slopes in 
these two groups, where cr is the standard deviation of the mean values. 



As mentioned in iNgeow & Kanburl ( 20I0l) . these six sets of P-L slopes can be 



grouped into two groups: one group with "shallow" slopes around -3.18, and another 
group with "steeper" slopes around -3.46. The expected IRAC band P-L slope can 
be calculated from = AnR^BxiT), hence Mjrac - -5 x gr \o%{P) + qt log(P) + 
const ant (where qr - 0.68 is the slope of period-radius relation taken from lGieren et al.l 



1 9991) ■ If assumed B;i{T) is constant at long wavelength, then the expected P-L slops is 
-5 X 0.68 = -3.40. On the other hand, due to Rayleigh-Jean approximation, B^{T) oc 
T at the IRAC band waveleng th, then aj ~ 0.18 estimated from color-temperature 
conversion (for details, see Neilson et al.l l2010). the expected P-L slope becomes -5 x 
0.68 -1-0.18 = -3.22. Both of these expected slopes are close to the averaged slopes in 
the observed "steep" and "shallow" groups. Finally, these empirical P-L slopes were 
compared to the slopes from synthetic P-L relations, based on theoretical pulsational 
models with varying metallicity, in Figure |2] Details of these synthetic P-L relations 
will be given elsewhere (Ngeow et al. - in preparation). The trends of synthetic P-L 
slopes with 12 -i- \og{0/H) > 7.9 appear to be in agreement to the empirical slopes 
given in "shallow" group. 
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Figure 2. Comparison of the six sets of empirical IRAC band P-L slopes to the 
synthetic P-L slopes from theoretical models (crosses). 



3. Ground-Based Light Curves for Cepheids Located Within Kepler Field-Of- 
View 

Kepler is a NASA space mission to search for Earth-like exo-planets using the transit 
method. Hence, the field-of-view for Kepler, with ~ 105deg^, is constantly pointing 
toward 19'^22™:-i-44'^30'. In addition to searching for exo-planets, light curve data from 
almost un-interrupted observations from Kepler are also very useful and valuable for 
asteroseismic and stellar variability study. Therefore, a consortium called Kepler As- 
teroseismic Science Consortium (KASC3 was formed. KASC consists of 14 working 
groups (WG), and WG7 is dedicated to Cepheids study with Kepler data. 

Since Cepheids are radially pulsating variable stars with periods longer than a day, 
it is not necessary to continuously monitor Cepheids in a given night or from different 
observatories across the globe. These multiple-site observing campaigns, on the other 
hand, have been practiced for other kinds of variable stars such as short period delta- 
Scuti stars to investigate, for example, various pulsating modes of these stars. Continu- 
ous photometry from Kepler, however, offers a gre at opportunity to study the possible 
non-radial pulsating modes as suggested by theory (Mulet-Mar quis et"al]|2007h . or even 
stochastically excited modes, in Cepheids. However, &;?/er observations are conducted 



jhttp : //astro . phys . au . dk/KASC7 
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in a customized broad-band filter, and complementary ground-based multi-color photo- 
metric and spectroscopic follow-up observations are needed for the Cepheid candidates 
located within the Kepler field-of-view. Detailed investigation of these candidates us- 
ing both of the Ke pler light curves and ground-based follow-up data can be found in 
Szabo et all (1201 ih . only a brief overview will be presented in this proceeding. 



3.1. BVRI FoUow-Up Observations for Selected Cepheid Candidates 

Prior to the launch of Kepler, VI 154 Cyg was the only Cepheid located in Kepler's 
field, while V2279 Cyg was a strong Cepheid candidate. Several other Cepheid candi- 
dates werealOT_selected based on the Kepler Input Catalog or other variable stars cata- 
logs ( Szabo et"aril201lh . Ground-based BVRI follow-up observations for these Cepheid 



candidates were carried out at Lulin Observatory (located at central Taiwan) and Tena- 
gra II Observatory (located in Arizona, USA), mostly from March to August 2010. The 
author is responsible for the observations from these two observatories, where the tele- 
scopes and CCD used are listed in Table 1. Figure[3]shows the distributions of FWHM 
in BVRI band images taken from these telescopes. Data reductions were done based on 
the following steps: 



Table 1 . Characteristics of telescopes and CCD for the BVRI follow-up observations. 



Tel. Abbreviation 


Observatory 


Aperture [m] 


CCD 


Pixel Scale ["/pix] 


LOT 


Lulin Obs. 


LO 


PI-1300B 


0.52 


SLT 


Lulin Obs. 


0.4 


Apogee U9000 


0.99 


TNG 


Tenagra II Obs. 


0.8 


STIe 


0.86 




Figure 3. Histograms of the FWHM distributions in BVRI bands, separated by 
the telescopes used. 
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1. Images were bias-subtracted, dark-subtracted and flat-fielded using IRAI^ Fringe 
corrections were applied to the / band images from LOT. 

2. Astrometric corrections were performed using astrometry . rtet ( Lang et al.l2010h 
and SCAMP (,Bertin.,20d6h . 

3. Instrumental magnitudes were extracted from the images using SExtractor 
(Bertin & Arnouts,1996) . These magnitudes were based on MAG_AUTO aperture 
photometry. 

4. Observations fr om Tenagra II O bservatory on 12 May 2010 included the Landolt 
standard fields ( Landolt 20091) . Instrumental magnitudes from this night were 
calibrated to the standard magnitudes. 

5. Photometry from other nights were relatively calibrated using the calibrated mag- 
nitudes from 12 May 2010. 

The resulted BVRI light curves for two of the candidates, VI 154 Cyg and V2279 Cyg, 
are present ed in Figure JH A dditional data were available from Sonoita Research Ob- 
servatorjH dSzabo et al.ll201 ih . 




Figure 4. Calibrated BVRI light curves for VI 154 Cyg (left panel) and V2279 
Cyg (right panel). Each light curves contain more than 70 data points. 

Light curves of VI 154 Cyg from Kepler o bservations resemb le a typical sawtooth 
shape of a Cepheid's light curve (see Figure 4 in lSzabo et al.ll201lh . Frequency analysis 



of almost continuous Kepler light curves for this variable showing a strong peak at the 
funda mental frequency, with other detectable harmonics in the spectrum (see Figure 
12 in lSzabo et al.ll201lh . This suggests that VI 154 Cyg pulsates radially in a regular 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Asso- 
ciation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 

"* Observations and the data reduction from Sonoita Research Observatory have been taken care by A. 
Henden. 
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fashion, without any non-radial or stochastic modes. The BVRI light curves of VI 154 
Cyg also strongly support the Cepheid nature of this variable. Fourier parameters (Rn 
and (pii), based on the V band light curve, of this Cepheid also fall within the distribu- 
tions defined by Galactic Cepheids. The phase lag from V band light curve and radial 
velocity curve confirms that VI 154 Cyg is a fundamental mode Cepheid. 

Even though the BVRI band light curves for V2279 Cyg mimic the light curve of 
a first overtone Cepheids, Fourier analysis of this variable revealed that this star is not 
a Cepheid. Furthermore, flares show up in Kepler light curves, and the light curve mor- 
phology is in agreement with rotational modulation. Additional spectroscopic follow- 
up observations confirmed that this variable is not a Cepheid. Cepheid nature of other 
candidates has been ruled out based on Kepler light curves and/or BVRI light curves, 
and V1154 Cyg remains the only Cepheid located within the Kepler's field. 

4. Conclusion 

I have presented some recent work on Cepheid research that carried out at NCU. In 
summary: 

• Mid-infrared P-L relations will be important in near future and in JWST era with 
the potential to reduce the systematic error of Hubble constant. Consequently, 
P-L relations based on Cepheids in Magellanic Clouds have been derived using 
the archival Spitzer data. 

• Almost continuous observation from Kepler is very valuable for stellar variability 
and pulsation study. However, ground-based follow-up observations are needed 
to complement Kepler data to fully investigate the variable stars in the Kepler's 
field. Together with Kepler light curves, the ground-based follow-up observations 
for Cepheid candidates confirmed that V1154 Cyg is the only Cepheid located 
within Kepler's field. 
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